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T T" 1 ' Abstract 
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-xt- : 

We discuss some aspects of recent research as well as more general issues about 

O ■ 

£f) • motivation, useful methods and open problems in the field of cosmological singu- 

| larities. In particular, we review some of the approaches to the general area and 

include discussions of the method of asymptotic splittings, singularity and com- 
pleteness theorems and the use of the Bel-Robinson energy to prove completeness 
^jrjl theorems and classify cosmological singularities. 

'Knowledge should be shared', yes, of course it should, but it is extremely 
doubtful whether an organic chemist would put missionary zeal into an at- 
tempt to interest the pure mathematician in his latest researches, or vice 
versa. To be realistic, the sharing of knowledge throughout the whole field 
of science is a pious hope, not a passion, and it is only passion that could 
overcome the inescapable difficulties of communication. J. L. Synge. 



1 Introduction 

Singularities have always been important, and all successful theories of physics have 
them by necessity. From the Newtonian dynamics of particles and fields, to relativistic 
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problems, to quantum dynamics, the ever-presence of singularities can for the most part 
be adequately treated one way or another, even avoided or ignored, although in some 
cases their consideration has triggered fundamental changes. 

But it is different in cosmology. Here, even though there are areas where the con- 
sideration of singularities can for the most part or wholly be neglected, no fundamental 
discussion of the basic cosmological issues can go on without facing them sooner or 
(preferably!) later. Up to the recent past, discussions of cosmological singularities have 
revolved around the same philosophical issues that we all know and grew up to respect, 
but which after a while any young, inquiring and productive cosmologist would wish to 
by-pass. But with certain exceptions. Not all discussions about singularities in cosmol- 
ogy are philosophical, they can be very mathematical indeed and at the same time touch 
upon very basic physical issues, probably at the cost of creating a situation that many 
people (even the creators themselves!) would wish to overpass. 

We believe that there lies ahead of us a golden era for research in cosmological 
singularities. The basic pillars are now almost done, the proliferation of newer and more 
interesting cosmological models is a call to arms, the mathematical techniques are better 
understood and explained and are definitely there waiting to be used, and cosmological 
observations, although richer than ever, do not constrain the possible models to any 
considerable degree. It is unfortunate that all this richness is often stated as a source for 
despair, it is a paradise for those who don't expect ready-made answers. 

In this review we are set to present the subject of cosmological singularities from 
a new angle, that of a logical and simplest form which however, may not be the most 
psychologically natural one. For motivation purposes, discussions of a new but already 
interesting subject are perhaps equally important to precise definitions and results, at 
least for the novice, and we have tried to present both here (resulting unfortunately in 
the usual unlucky compromise). In the next section after a short review of the basic 
definitions of completeness and singularity, we introduce the problem of cosmological 
singularities in two stages: Firstly, it is set at the level of special, exact cosmologies that is 
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when we have chosen a particular spacetime, a specific matter component and a geometric 
action to describe the gravitational field, for instance a fluid-filled Friedmann universe 
in general relativity. Secondly, we formulate the singularity issue at the different level of 
leaving the choice of spacetime unspecified. This leads us in turn into two ways of attack, 
one through spacetime differential geometry to formulate and prove general results about 
geodesic incompleteness (singularity theorems), the other being the analysis of the field 
partial differential equations describing the evolution of Einsteinian spacetimes. 

It is important to realize that any information about the singularities of a system in 
the general case (that is without assuming any symmetry) cannot be usefully reduced to 
give us information about particular instances of the field in situations of symmetry (and 
vice versa of course!). This last case requires a different analysis. A concrete general 
technique, the method of asymptotic splittings, designed specifically to deal with the 
approach to a spacetime singularity of exact cosmologies is presented in section 3. The 
aim of this method is to decompose the vector field defining the model in a suitable 
way and follow the dominant features of it all the way down the finite time singularity. 
This leads to the construction of asymptotic series expansions of the unknowns in the 
neighborhood of the singularity, from which valuable information about the field can be 
extracted. Some examples of the application of the method of asymptotic splittings are 
discussed. This method is an adaptation to real dynamical systems (finite dimensions) 
of the philosophy of certain blow-up techniques in partial differential equations, and of 
methods that exist in the field of complex differential equations (singularity analysis). 

Next we move on, in section 4, to discuss sufficient and necessary conditions for 
the occurrence of singularities. This section, besides the classic singularity theorems, 
includes more recent results about singularities associated primarily with the name of Y. 
Choquet-Bruhat. We also provide applications of these criteria for singularity formation 
to the isotropic category and show how these lead to a novel classification of spacetime 
singularity types therein through the use of the Bel-Robinson energy. We also use the 
Bel-Robinson energy to give a new characterization of the g-completeness and eternal 
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acceleration of Robertson- Walker (RW) spacetimes. 

The absence of a discussion section at the end is a healthy one, and conclusions are 
out of place at this initial stage of development of this most beautiful subject. 

I would like to thank many colleagues and collaborators, especially Ignatios An- 
toniadis, John Barrow, Yvonne Choquet-Bruhat, Ifigeneia Klaoudatou, John Miritzis, 
Vincent Moncrief, Alan Rendall and Antonios Tsokaros, whose work and ideas I have 
freely used in this review, for many stimulating discussions about various important 
issues associated with spacetime singularities. Part of the research reviewed here was 
conducted while the author visited CERN-Theory Division in Geneva and IHES in Bures 
sur Yvette and I thank those institutions for allowing me to use their excellent facilities. 

2 Approaches to the singularity problem 

In this section we give the basic definitions of what it means for a spacetime to be 
singular. We then discuss the general problem of spacetime singularities in the framework 
of exact cosmological solutions (physical approach). Finally, we extend the formulation 
of the problem to the case of spacetimes without symmetries and discuss various more 
mathematical approaches used to tackle this most basic cosmological issue. 

2.1 Definitions 

We give in this subsection the basic definitions of complete as well as singular space- 
times. Let (M,g) be a spacetime with M. = (a, b) x E, where (a, b) is an interval in 1R 
(possibly the whole line), and the spatial slices E t (= E x {£}) are spacelike submanifolds 
endowed with the time-dependent spatial metric g t . The standard geometric definition 
of spacetime completeness in general relativity is the following. 

Definition 2.1 (Geometric) We say that the spacetime (Ai,g) is non singular if it 
is causally geodesically complete, that is every causal (i.e., timelike or null) geodesic 
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7 : (a, b) C R — > 7W , defined in a finite subinterval, can be extended to the whole real 
line. 

However, this definition is not useful in evolution problems because in such studies one 
is typically interested not in the behaviour of geodesies but in that of the dynamical 
variables of the problem. These satisfy differential equations such as the Einstein or 
some other field equations. In such problems one is typically interested in proving local 
or global in time existence of the spacetime in question. This leads us to the following 
alternative. 

Definition 2.2 (Dynamical) We say that the spacetime (M.,g) exists globally as a 
solution of the field equations if 

1 1 (g t , dgt, 0) | \lp < oo, uniformly for all p, t. 

This definition, deliberatively vague, means intuitively that a spacetime such that its 
spatial metric variables g t , their derivatives dg t and any matter fields present are finite 
for all time, must be considered as globally existent. In a specific problem of course one 
must choose the appropriate norms so that they make sense with the given conditions on 
the various data. However, even when, as it is usually the case, global existence comes 
together with non singularity (in the above sense of geodesic completeness), one must 
prove that this is indeed the case as the two notions are different. 

The negations of the above definitions lead to the notions of singular and finitely 
existent spacetimes. 

Definition 2.3 (Geometric) We say that (A4,g) is singular if it is geodesically in- 
complete, that is there is at least one incomplete causal geodesic 7 : (a, b) C R — > A4. 

For a spacetime to have only a finite duration, we require the following 

Definition 2.4 (Dynamical) We say that (A4,g) exists as a solution of the field equa- 
tions for only a finite time if the following holds: 

1 1 (gi, <9g t , 0) I \lp becomes unbounded at some p, t. 
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These definitions allow us to make several remarks. 

Remark 2.1 (Equivalence) The two definitions, geometric and dynamical, of com- 
pleteness on the one hand, and singularity on the other, are considered as equivalent, but 
there is to our knowledge no formal proof of this in the literature. So one must check 
for instance that a spacetime that exists globally for all future time is also geodesically 
complete to prove that it contains no singularities. For an example of this, see |3J Q]/. 

Remark 2.2 (True singularity) It is important to emphasize that one should be care- 
ful not to mix a singularity in the congruence of a bunch of geodesies with a true spacetime 
singularity. The singularity theorems (see below) give conditions when one ends up with 
a true spacetime singularity, the proof, however, of such statements is by contradiction 
and therefore one cannot use such arguments to really construct incomplete spacetimes. 

Remark 2.3 (Fixed and sudden singularities) The actual finite duration of a sin- 
gular spacetime (according to the dynamical definition), in other words the position on 
the time axis of the finite time singularity, is in fact its maximum interval of existence 
considered as a solution. This therefore may very well depend on the initial and final 
conditions set in the problem. One is therefore led to consider fixed singularities, those 
which have positions independent of the initial conditions as well as movable, or sudden 
singularities, whose position is actually very sensitive to the initial conditions. An in- 
teresting nontrivial problem is to characterize these two types of singularities using the 
geometric definitions of geodesic incompleteness. For more details see Section \Ji~JZ . 

2.2 The physical approach 

Up to around the 1960's only three dynamical theories existed in physics: Newtonian 
dynamics (of particles, continua and fields), relativistic dynamics (including both the 
special and the general theory) and quantum dynamics (splitting naturally in its new- 
tonian (mechanics) and relativistic (field theory) parts). Since then, however, there has 
been a proliferation of different dynamical theories some aiming at rectifying one or more 
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of the issues which arose in one of the three pillars, others pointing at the other direction, 
that of the unification (eventually of relativistic and quantum dynamics). This sparkling 
diversity is especially apparent in theories which include gravity, and the following table 
aims at offering a broad but humble picture of this most chaotic situation. 



Cosmologies 



Theories of gravity 


Spacetimes 


Matterfields 


General Relativity 


(Anti-) de Sitter 


Vacuum 


Higher Derivative 


RW 


Fluids 


Scalar- Tensor and VC 


Bianchi, 


Scalar fields 


Inflation and QC 


GSdel 


n-form fields 


Braneworlds 


T-NUT-M 


Phantoms, tachyons 



With the advent of string theory and its subsequent generalizations and extensions pre- 
vious dynamical theories of gravity found their natural places under its umbrella in the 
form of effective actions and field equations and their overall position gave to the picture 
a notable matter-of-fact-ness. 

But it is disturbing. The whole picture is blurred with the existence and apparent 
(mathematical at least) viability of all these theories, while their structures are connected 
to one another by complicated laws, duality and conformal transformations, making it 
impossible to choose some and to leave more out (the existence of the so-called string and 
Einstein frame formulations are good examples of this). We don't yet know what really 
lies behind this web of theories but the real question is: How is it that we can tolerate 
this apparently arbitrary multiplicity? The standard answer to this pressing question is 
that of course each and every one of these possible theories in the web represent different 
viable mathematical models of reality. 

The least thing one needs in this game is to become inflated with long words. The 
singularity problem acquires the following simple form in this context: Analyze the pos- 
sible models (like those that can be readily extracted from the Table above) and discover 
what they have to say about the existence and nature of their spacetime singularities. It 
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is to be expected that the obvious proliferation of models will lead to an equally obvious 
number of different kinds of singularities. This is a good thing, it is like discovering new 
kinds of life in zoology (for more on this, see [HE])- 

The essential mathematical tools needed to construct such cosmological models (with 
the aim eventually to treat their singularities) come from tensor analysis. One ends up 
with tensor equations, invariant under arbitrary coordinate transformations, and these 
describe the dynamical content of the model. Then, techniques from the theory of dy- 
namical systems enter the picture, and with geometry and dynamical systems one hopes 
to be able to tackle the many different subtleties which will be pertinent to the singularity 
structure of the spacetimes in question. Due to its richness, unraveling the singularities 
of different cosmological models is a game that must be played with adroitness, many 
surprises are still hidden from our ever expanding horizon. More on this approach will 
be given in Section 3. 

2.3 The mathematical approach 

There is another way of looking at the problem of cosmological singularities, but before 
we discuss this in this subsection, let us summarize in a sentence the content of the 
approach taken in the previous subsection. For the admitted variety of theories of gravity 
we have to pay a price, namely, the spacetimes to-be-used are mathematically simple, 
perhaps too simple, but there are at least two good reasons that make that a legitimate 
approach: Even if one does not believe that the actual models constructed this way 
are true in reality, one still has a reason to consider them for only in this way we have 
even the slightest chance to use the available observational data (cf. [6], p. 470). Or, 
we perhaps, like Einstein, must elevate the homogeneity assumption to the status of a 
physical principle, cf. [7], p. 16). Follow either of these two leads and you conveniently 
find yourself walking along with the other physical cosmologists. 

Others, however, follow none. The point of view of the mathematical cosmologist, 
at least in its traditional sense, is that, because we know so little from observations, it 
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is legitimate to consider arbitrary manifolds and metrics (cf. [8], p. 30). There is also 
another reason for this legitimacy. It is quite possible that homogeneities at small scales 
might change drastically the global mode of behaviour in any model when we follow its 
evolution in either early or late (proper) times. For the treatment of homogeneous cos- 
mological models in any theory of gravity, this could affect the results obtained through 
the physical approach quite seriously, even deadly^. 

The mathematical cosmologist does not have to go far to find problems of his liking, 
general relativity is for the most part adequate. Even the simplest of theories when looked 
from such a perspective can have a richness of unprecedented beauty. However, sooner 
or later the mathematical approach faces a true obstacle: The transcendental difficulty 
of the equations. The Einstein equations are, in their most interesting form, a system 
of hyperbolic evolution equations together with elliptic constraints and this system will 
obviously have to be augmented by further equations describing the evolution of the 
various matter fields. 

However, some progress has been made and here we have the Cauchy problem (see [8] 
and references therein). In this approach, spacetimes are not known in advance, they are 
in fact built up from initial data as the solutions that describe them evolve, or develop, 
through the field equations from that data. The issues of the existence and uniqueness 
of maximal developments, as well as that of the nature of singularities in these solutions 
are the central objects of study here, and it is of interest to apply these methods not 
only to general relativity but also to other interesting theories of gravity formulated in 
the string frame, something that has not yet been tried to any measurable degree of 
completeness. 

Another area in which modern mathematical methods can be of great interest is 



1 There are of course certain 'no-go' conjectures, the BKL conjecture being the primary example (see, 
the very recent work [9] for an account of the current status of the proof of this conjecture). For the 
singularity problem, this basically says that what happens in the most general homogeneous model on 
approach to the singularity (technically for the Bianchi IX spacetimes) describes the generic situation, 
that is the approach to the singularity when there are no symmetries. 
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in fact that which contains models of interest to the physical cosmologist. Here a more 
mathematical approach to models already formulated by physicists can help in improving 
the rigor of many presentations as well as unify and polish the essential elements of many 
models considered in the literature. We shall have to say more on this later. 

Still another direction to the singularity problem open to the modern mathematical 
cosmologist is through Lorentzian geometry. This approach has led historically to the 
first singularity theorems, those of R. Penrose and S. W. Hawking (see [TDJ for a modern 
introductory — of these thecal Recently there has been a related body of 
results, the completeness theorems (see jl] and also [12] for an application of these results 
to cosmology). These theorems aim to give sufficient conditions for a spacetime to be 
geodesically complete, however, like the singularity theorems (but unlike the dynamical 
mathematical approaches described above), they require a spacetime of a known form in 
order to be applied. We shall have to say more about these theorems in Section 4. 



3 Geometric asymptotics for exact models 

A large part of the literature^ on cosmological singularities is directly concerned with 
the problem of tracing out the exact, early or late time, asymptotic behaviour of various 
cosmological models. In all these studies, a cosmology (in the sense of the previous 
Section) is given and the problem is then to determine the nature of the spacetime 
and of the various fields in the neighborhood of a past or future singularity. This is 
basically an asymptotic problem since the form of the spacetime is known in advance 
but the precise asymptotic relations of the unknown quantities present in the cosmology 



2 There are many results in this and related directions, see the book [IT] where this whole approach 

is presented along with interesting theorems and open problems. 

3 To do full justice to the enormous number of published papers in the many different areas of this 

subject can lead to a project in itself, something that this humble review never intended to do (a 

comprehensive Bibliography of Papers on Cosmological Singularities is currently under preparation). 

We apologize in advance to many readers to whom the choice of references included here, and more 

importantly their omissions, seems arbitrary. The present author is aware of this issue only too well! 
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are to be determined. In this Section, we present an outline of a new (or more precisely 
cleaned-up) method, called hereafter the method of asymptotic splittings (cf. |T3j for a 
more detailed presentation), which enables us to precisely determined such asymptotic 
relations. The utility of the method of asymptotic splittings to analyze cosmological 
singularities is then shown in a number of recently examined cases (with references to 
the relevant papers). We also add some comments and remarks about the issue of 
fixed vs. spontaneous (or sudden) singularities, the latter having recently gained some 
attention, in a slightly different context, through the work of J. D. Barrow and others 
(see [TH [151 [12] an d refs. therein). 

3.1 Method of asymptotic splittings 

In this subsection, we follow (13J closely and this reference is to be consulted for notation, 
explanations, and proofs of any stated result. The general setting is a vector field f : 
M. n — > TA4 n and the associated dynamical system defined by f on the manifold Ai n , 

x = f(x), (3.1) 

with (■) = d/dt. Depending on the number of arbitrary constants, a solution of our dy- 
namical system can be general, particular or exact. Any solution, however, can develop 
finite-time singularities, that is instances where a solution x(i; c\, ■ ■ ■ , C&), k < n, mis- 
behaves at a finite value t* of the time t. This is made more precise as follows. We say 
that the system x = f (x) (equivalently, the vector field f) has a finite-time singularity if 
there exists at t G K and a x G Ai n such that for all M e R there exists an 5 > such 
that 

||x(t;xo)||i* > M, (3.2) 

for \t — i*| < 8. Here x : (0,6) — ► Ai n , x = x(i ) for some t Q G (0,6), and || • \ \lp is 
any L p norm. We say that the vector field has a future (resp. past) singularity if > t 
(resp. t* < to)- Note also, that to is an arbitrary point in the domain (0, 6) and may be 
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taken to mean 'now'. We often write 

lim ||x(*;x q )||lp = oo, (3.3) 

to denote a finite-time singularity at £*. Our basic problem then is to find the structure 
of the set of points x in A^ n such that, when evolved through the dynamical system 
defined by the vector field, the integral curve of f passing through a point in that set 
satisfies property (13.31) . 

There are two kinds of finite time singularities that a nonlinear dynamical system can 
possess, fixed and movable (sudden, or spontaneous in other terminology is often used 
for a movable singularity). A singularity is fixed if it is a singularity of x(£; C) for all C; 
otherwise, we say it is a movable singularity. Our basic problem is then: What can a 
vector field do, or equivalently, how do the solutions of the associated dynamical system 
behave in the neighborhood of a finite time singularity? Assume that we are given a 
vector field and suppose that at some point, t*, a system of integral curves, corresponding 
to a particular or a general solution, has a (future or past) finite-time singularity in 
the sense of definition H 3 . 2 f) . The vector field (or its integral curves, solutions of the 
dynamical system defined by C) can basically do two things sufficiently close to the 
finite-time singularity, namely, it can either show some dominant feature or not. In 
the latter case, the integral curves can 'spiral' in some way around the singularity ad 
infinitum so that (13.21) is satisfied and the dynamics is totally controlled by subdominant 
terms, whereas in the former case solutions share a distinctly dominant behaviour on 
approach to the singularity at £* determined by the most nonlinear terms. 

Our approach to this problem is an asymptotic one. We decompose, or split, the vec- 
tor field into simpler, component vector fields and examine whether the most nonlinear 
one of these shows a dominant behaviour while the rest become subdominant in some 
exact sense. We then built a system of integral curves corresponding, where feasible, to 
the general solution and sharing exactly its characteristics in a sufficiently small neigh- 
borhood of the finite-time singularity. In this way we are led to a general procedure to 
uncover the nature of singularities by constructing series expansion representations of 
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particular or general solutions of dynamical systems in suitable neighborhoods of their 
finite-time singularities. 

This method of asymptotic splittings [TH] consists of building splittings of vector fields 
that are valid asymptotically and trace the dominant behaviour of the vector field near 
the singularity. A resulting series expansion connected to a particular dominant balance 
helps to decide whether or not the arrived solution is a general one and to spot the exact 
positions of the arbitrary constants as well as their role in deciding about the nature of 
the time singularity. To apply the method of asymptotic splittings to a given dynamical 
system so as to discover the nature of its solutions near the time singularities, we must 
follow a number of steps: 

1. Write the system of equations in the form of a dynamical system x = f(x) with 
x = (xi, • ■ ■ , x n ), and identify the vector field f (x) = (/i(x), ■ ■ • , / n (x)). 

2. Find all the different weight-homogeneous decompositions of the system, that is 
the splittings of the form 

f = f(o) +f (D + ... + f (fc) ) 

and choose one of these splittings to start the procedure. 

3. Substitute the scale-invariant solution 

x (°)(r) = ar p , 

into the equation x = f^. Study the resulting algebraic systems, and find all 
dominant balances (a, p) together with their orders. 

4. Identify the non-dominant exponents, that is the positive numbers 
such that 

fub,Q-) (rP) 

T H ~ ► U. 

r p-l 
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5. Construct the K-matrix /C: 

f(°) _> £>f(°) ^ J Df(°)( a ) -> Df (0) (a) - diagp. 

6. Compute the spectrum of /C, 

spec(/C) = (-l,p 2 , • ■ • ,p n )- 
Is /C semi-simple? Are the balances hyperbolic? 

7. Find the eigenvectors of /C. 

8. Identify s as the multiplicative inverse of the least common multiple of all the 
subdominant exponents and positive K-exponents. 

9. Substitute the Puiseux series 

into the original system. 

10. Identify the polynomials Pj and solve for the final recursion relations which give 
the unknown coefficients Cj. 

11. Check the compatibility conditions at the K-exponents, 

vj • P p = 0, for each eigenvalue p. 

12. If the Puiseux series is valid, then the method is concluded for this particular 
splitting. Otherwise, if compatibility conditions are violated at the eigenvalue p*, 
restart from step M by substituting the logarithmic series 

too oo 
a + Yl Ci ^ i/S ( rP log r ) j/S 
i=l j=l 



oo 

i=o 
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13. Get coefficient at order p*. Write down the final expansion with terms up to order 



14. Verify that compatibility at p* is now satisfied. 

15. Repeat whole procedure for each of the other possible decompositions. 

Following the above steps even up to that of calculating a dominant balance in one 
particular decomposition (that is up to step 4), can be very useful since it offers one 
particular possible asymptotic behaviour of the system near the time singularity. In this 
respect, the whole method expounded here is truly generic since it helps to decide the 
generality of any behaviour found in an exact solution - that is, how many arbitrary 
constants there are in the final solution that shares that behaviour (particular or general 
solution). It is rare that a Puiseux series is inadequate to describe the dynamics (semi- 
simplicity of K), but in such cases one must resort to the more complex logarithmic 
solutions. 

3.2 Example: Gauss-Bonnet cosmologies 

There are many interesting systems that could be seen under the magnifying glass of 
the method of asymptotic splittings, cf. the systems in [131 UHl HH1 [20] . Here we present 
a short summary of some of the results of our very recent work [21] in collaboration 
with A. Tsokaros, which is to be followed for motivation and further explanations. The 
starting point is the action 



Jm 4 

where £ to tai * s the lagrangian density of the general quadratic gravity theory given in the 

forilj] £ t otal =£(R) + ^matter, with 





(3.4) 



C(R) = R + BR 2 + CRic 2 + DRiem 2 , 



(3.5) 



the conventions for the metric and the Riemann tensor are those of |17j . 
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where B,C,D are constants. Since in four dimensions we have the Gauss-Bonnet identity, 

6 / R 2 GB dfi g = 0, R 2 GB = R 2 - 4Ric 2 + Riem 2 , (3.6) 

in the derivation of the field equations through variation of the action associated with 
( 13. 5p . only terms up to Ric 2 will matter. Below we focus exclusively in spatially flat 
universes of the form 

ds 2 = dt 2 - b{t) 2 {dx 2 + dy 2 + dz 2 ), (3.7) 

which are radiation dominated (P = p/3), and use only the 00-component of the field 
equations, this being the following equation satisfied by the scale factor b(t): 

b 2 



b 2 K 



2 i^ + 9^---3^ 

b 2 b 3 b 2 ¥ 



b 2 



where b\ is a constant defined by 



2 



8nGp b\ 



(from ViT i0 = 0). (3.9) 



3c 4 6 4 ' 

Note that the Friedmann solution \/2b\t of general relativity satisfies the above equation. 

To apply the method of asymptotic splittings, we set b = x, b = y and b = z, and 
then Eq. (I3.8P can be written as a dynamical system of the following form: x = f(x), 

x = (x,y,z): 



x = y, y = z, z = — - 7 —- 2 - — + — + — (3.10) 



y_ _ b'j _V± + ^_ + ^_ 
2k 2nyx 2 x 2y 2x 2 
Then we look for the possible weight-homogeneous decompositions. The most interesting 

one has dominant part 

fV = (y,z,^-^ + 3 £), (3.ii) 



2y x 2x 2 



while the subdominant part reads 



f! " b =(°'°'-^+l;)' (312) 

with f = f(°^ + f sub . The dominant balance for this decomposition turns out to be the 
unique (full order) form 

cp> =(64 -i). ' 313 > 
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where a is an arbitrary constant. 

Next we turn to an examination of the form and properties of the K-matrix. The 
Kowalevskaya exponents for this particular decomposition, eigenvalues of the matrix /C = 
Df (a) — diag(p), are {—1, 0, 3/2} with corresponding eigenvectors {(4, —2, 3), (4, 2,-1), 
(1, 2, 2)}. (The arbitrariness of the coefficient a in the dominant balance reflects the fact 
that one of the dominant exponents is zero with multiplicity one.) After following the 
other steps of the method we finally arrive at the following series expansion representation 
of the behaviour of the system in the neighborhood of the initial collapse singularity: 

x(t) = a (t - to) 1 / 2 + c 13 (t - t ) 2 + ^f^r (* - *°) 5/2 + • ■ ■ • ( 3 - 14 ) 

The series expansions for y(t) and z(t) are given by the first and second time derivatives 
of the above expressions respectively. 

There are several comments to be made for this solution. First, the series (I3.14p has 
three arbitrary constants, a, C13, to (the last corresponding to the arbitrary position of the 
singularity) and is therefore a local expansion of the general solution around the movable 
singularity t . Secondly, our solutions are stable in the neighborhood of the singularity 
in the sense that all flat, radiation solutions of the quadratic gravity theory considered 
here are Friedmann-like regardless of the sign of the R 2 coefficient. This implies that 
there are no bounce solutions of the model analyzed here. It remains an open problem to 
find a physically reasonable modification of this model with stable, bouncing solutions 
valid at early times. 

4 Differential geometric issues 

There are many instances where one does not wish to put precise assumptions on the form 
of the cosmological spacetime but leave it unspecified instead. In such circumstances it 
is important to be able to find general conditions and criteria under which a 'universe' 
represented in such a generic way will develop a physical, past or future, singularity in a 
finite time. In fact, it turns out that one can formulate precise criteria for such purposes 
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and even begin to classify the possible mathematical and physical singularities of this 
generic situation. 

In this section we review some attempts to tackle this problem, beginning with the 
(statement of the) time-honored singularity theorems, results predicting the general ex- 
istence of spacetime singularities, first formulated in the sixties. These results utilized in 
their proof some of the most beautiful concepts and techniques from differential topology 
and Riemannian geometry, but since then the general theory has remained somewhat 
dormant and perhaps wanting. Recently, it has been possible to revive interest in this 
whole area and new results predicting the existence of complete as well as singular 
spacetimes have been made possible, mainly due to the efforts of Y. Choquet-Bruhat. 
An adaptation of these ideas to simple cosmological solutions reveals structures which 
are perhaps difficult to unearth by more traditional methods and analysis. 

4.1 Completeness and singularities 

The geometric criteria to be developed in this subsection for a spacetime to be either 
geodesically complete or incomplete will be based on certain properties of its spacelike 
submanifolds. Below we focus exclusively on spacelike hyper surf aces, although subman- 
ifolds of higher codimension also play a significant role for such purposes (the most 
well-known example of the latter situation is the notion of a closed trapped surface). 

Consider a spacetime of the form (Ai,g) with M. = £ x X, X being an interval in R 
(for simplicity we can take it to be the whole real line) and (£,<?*) a smooth, spacelike 
submanifold of dimension n, in which the smooth, (n+l)-dimensional, Lorentzian metric 
g splits as follows: 

g = -N\9°) 2 + g t , g t = 9ij Q l Q\ 9° = dt, 9 l = dx l + f3 l dt. (4.1) 

Here N = N(t,x l ) denotes the lapse function, (5 l {t,x^) is the shift vector field and all 
spatial slices S( (= Sx {t}) are spacelike hypersurfaces endowed with the time-dependent 
spatial metric g t = gijdx l dxK We assume that our spacetime (V, g) is regularly sliced, 
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that is the lapse function is bounded, the shift vector field is uniformly bounded and the 
spatial metric is itself uniformly bounded below, cf . pQ [2] . 

There are two kinds of information about the slices £ t that play a role in what follows: 
topological and geometric. The topological part of the theory to-be-developed starts with 
a connection of the completeness of Y, t with the global hyperbolicity of (Ai, g), as in the 
following result from [TJ, [2] which may be consulted for the proof. 

Theorem 4.1 Let (Ai,g) be a regularly sliced spacetime. Then the following are equiv- 
alent: 

1. (So, go) is a go-complete Riemannian manifold 

2. The spacetime (A4,g) is globally hyperbolic 

We thus see that our ability to successfully exercise classical determinism in a general 
spacetime is intimately connected to the topological property of geodesic completeness, 
not of the whole of spacetime, but of each one of its spacelike slices (considered as 
Riemannian manifolds, so that, in particular the Hopf-Rinow theorem is valid). Unfor- 
tunately to say something about the geodesic completeness of the whole of our general 
spacetime (A4,g) is not possible by using solely topological pieces of information as the 
one given above (unless (A4,g) is trivial in some sense, see[2]). 

We need to know in addition more about the geometry of (M., g). What we specifically 
need to know in order to state the following singularity as well as completeness theorems, 
is some submanifold geometry (cf., for example, [10] chap. 4 and parts of chap. 10). For 
both kinds of theorem, the notion the shape tensor (or second fundamental form) of the 
slice Et becomes vital. Consider the mean curvature vector field H of Tit, denote by V 
the induced connection on S 4 , and let K = norV be its shape tensor which is defined 
on pairs of vector fields of S t . 

The trace and norm of K play a fundamental role. For a unit future pointing vector 
field U normal to S t , we define the convergence of T, t to be the real- valued function 9 on 
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the normal bundle NTi t 

= (U, H) = — *— trace K. (4.2) 
n — 1 

A natural interpretation of 9 is to say that when 9 is positive the spacelike slice is bent 
inward so that its outward pointing normals converge, whereas when it is negative we 
have the reverse situation. One may easily then imagine a spatial slice such that 9 is 
positive everywhere on it: we will picture it as a convex spacelike arc of infinite extent. 

The only other piece of geometric information we shall need is this: the gt norms 
of K as well as of the spatial gradient of the lapse function, namely, \K\ St and |ViV| gt . 
These are the natural objects to consider and estimate in an effort to see how large or 
small they can become and so to develop a feeling on how the slice itself evolves with 
proper time. 

We are now ready to state a basic singularity theorem and a completeness theorem 
to get a taste of this kind of result (for more on this see the references above). The first 
result is Hawking's singularity theorem. 

Theorem 4.2 If: 

• Ric(X,X) > for all causal vector fields X of the spacetime (Ai,g) 

• 9 > C > 0, everywhere on the Cauchy slice S 7 

then no future- directed causal curve from £ can have length greater than 1/C. 

This theorem continues to hold if in the place of the words 'Cauchy slice' one substitutes 
the weaker 'compact slice', although in this case the proof is somewhat more involved 
(cf. [10J pp. 431-3). The proof is by contradiction: One cannot have both a complete 
spacetime and a Cauchy (or compact) slice with a positive convergence everywhere. 
That is unlike the case where the convergence were positive only locally, then nothing 
bad would happen. (This is the case, for instance, when one considers a convex spacelike 
arc of finite length with positive convergence in two-dimensional Minkowski space: its 
normals are converging, they meet at some point to the future of the arc but there is 
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no singularity, spacetime is geodesically complete.) Therefore this theorem predicts that 
under this situation spacetime must be geodesically incomplete, that is singular, but in 
a very special way: the whole of spacetime will eventually fold around into itself! 
On the brighter side, we have the following completeness result 

Theorem 4.3 If 

• (-M,g) is a globally hyperbolic, regularly sliced spacetime 

• for each finite ti, |ViV| gt and \K\ St are integrable functions on [ti,+oo) ; 
then (Ai,g) is future causally geodesically complete. 

The proof of this result is achieved by showing that all future-directed geodesies have 
an infinite length, and it is done by analyzing properties of the geodesic equation. For 
various applications of this theorem to cosmological spacetimes, see [T2l ITS] . 

4.2 Fixed and sudden singularities 

We pause in this subsection to discuss, in passing, some qualitative differences between 
the two possible dynamical kinds of finite time singularities inherent in the method 
of asymptotic splittings, namely, fixed and movable (or sudden) singularities. Both 
types of singularities are of course in accordance with the assumptions of the singularity 
theorems, although in discussions of the latter the distinction is never made explicit. 
(For concreteness, the reader if he so wishes, may take a simple cosmological metric, for 
instance one from the Robert son- Walker family, which has only one unknown function, 
the scale factor a.) We assume that we have a situation where the assumptions of 
the singularity theorems are valid and the model has a true, finite time singularity. Our 
problem currently is to obtain conditions under which this singularity is fixed or movable. 

Suppose we have a singly infinite family of geodesies in our spacetime. These geodesies 
form a 2-dimensional geodesic surface parametrized by 




(4.3) 
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and one understands that the geodesies defined in this way are the parametric lines of u 
in the geodesic surface. We denote by p r = dx r /du the unit tangent vector field to the 
geodesies, and by rf = (dx r /dv)dv the connecting, infinitesimal vector field between two 
adjacent geodesies of the family, also called a Jacobi field. These two vectors p r , rj r satisfy 
two fundamental equations, namely the geodesic equation and the Jacobi equation (or 
equation of geodesic deviation) respectively Both equations are very familiar and play 
a basic role in any discussion of geodesic incompleteness. They read: 

f + r mnP>™ = Geodesic Equation, (4.4) 

and 

V r + R r smnP s V m P n = 0, Jacobi Equation. (4.5) 

Here, differentiation is meant with respect to arc length (wherever it makes sense). 
Both equations are nonlinear but their main difference is this: Whereas the equation 
of geodesies is not a tensor equation but a frame-dependent one and so any conclusions 
using it may change under coordinate transformations, the Jacobi equation is a tensorial 
equation and anything we say using it will be valid in an invariant geometric way. 

For instance, although the existence of conjugate points (two events each having 
rj = 0) on two adjacent geodesies follows from an analysis of the Jacobi equation and is 
true in any system once it is established, properties of solutions of the geodesic equation 
are not so insensitive. 

Suppose we have a Robertson- Walker metric with a true singularity, as, for instance, 
the initial big bang (collapse) singularity in the standard cosmological models, cf. [17], p. 
366. The actual spacetime point, location of the singularity, is then a conjugate point in 
the sense that there is a non-zero Jacobi field between any later point and the big bang 
that becomes zero there. Timelike geodesies that end at the big bang (which in this 
case is an all-encompassing singularity) are still solutions of the geodesic equation (I4.4p . 
From standard expressions, the Christoffel symbols diverge at the singularity, in other 
words the coefficients in the geodesic equation ( 14.4ft are diverging. This means that the 
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standard big bang singularity is a fixed singularity in the terminology of the previous 
Section. It does not depend on the initial conditions (which may be taken to be at any 
later time), in this case the singularities of the dynamical system are read off from the 
singularities of its coefficients. 

A completely different situation may however arise, and this is reminiscent with 
what happens to solutions of the very simple nonlinear equation x = x 2 with initial 
condition x(0) = xq. There is the solution x(t) = x$/ (1 — x$t) and this is defined only on 
(—00, l/x ), when x > 0, on the whole of R for x Q = 0, and on (l/x , 00) for x < 0. We 
see that the solution ceases to exist after time 1/xq. The coefficient in the equation is, 
however, regular (in fact equals one). So this equation has singularities whose position 
depends on the initial conditions and they 'move' accordingly with the initial condition, 
but the coefficients of the equation don't 'see' them. This is the simplest example of a 
movable (or sudden) singularity. 

We believe that in general relativity we can in principle have a similar situation. 
For, if a geodesic defined on an interval (a, b) is, say, future incomplete (as that would 
follow from a tensorial analysis of the invariant Jacobi equation), then it cannot be 
extended after b to the whole of the real line, but this b may or may not depend on 
the initial conditions chosen, making the singularity either fixed or movable. This is 
only to be expected because the geodesic equation f)4.4p is nonlinear and so cannot have 
singularities of only the fixed kind (that is those that are read off from the coefficients 
of the equation). Singularities in the traditional form of g-incompleteness are usually 
thought of in the literature only as fixed singularities, no dependence of the singularity 
on the 'initial' (later for a past, earlier for a future movable singularity) conditions. 

4.3 Bel-Robinson energy 

It is clear from what has been said earlier that there are several different, inequiva- 
lent ways to characterize cosmo logical singularities. First, there is the traditional way, 
through differential-geometric methods, as explained in the standard text on this ap- 
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proach [22] . This set of methods uses a global approach to spacetimes with an emphasis 
on the causal, topological aspects of singularities. But there is more to characterizing 
singularities than this, and using dynamics one may arrive at a complementary set of 
tools for such purposes. We saw how, using the method of asymptotic splittings, we 
can have quantitative information about a found singularity and form a clear picture 
on the different modes of approach through series expansions of the key unknown func- 
tions pertinent to the given problem. Other relevant approaches based on dynamical 
systems are also of great use, cf. [23]. A third approach is to use the contrapositive 
of the completeness theorem discussed earlier and end up with necessary conditions for 
singularities. This approach to the nature of cosmological singularities was taken up in 
[T2l [T6] for the simple case of Robert son- Walker cosmologies. The result is that we can 
arrive at an interesting trichotomy of the singularities in various Friedman cosmologies 
by using the two main functions of the problem, namely, the Hubble expansion rate H 
and the scale factor a. 

To the three main approaches to singularities another must be added. In none of 
these does the presence of matter fields play a role, at least directly, and the following 
way to classify cosmological singularities uses exactly this missing piece of information: 
It is based on the introduction of an invariant geometric quantity, the Bel-Robinson 
energy which takes into account precisely those features of the problem, related to the 
matter contribution, in which models still differ near the time singularity while having 
similar behaviours of a and H. In this way, we arrive at a complete classification of the 
possible cosmological singularities in the isotropic case (for more details see [MIEHl 26J). 

The Bel-Robinson energy at time t is given by 

B(t) = ~ [ (\E\ 2 + \D\ 2 + \B\ 2 + |iJ| 2 ) d»- gt , (4.6) 

where by \X\ 2 = g^ g kl XikXji we denote the spatial norm of the 2-covariant tensor X, and 
the time-dependent space electric and magnetic tensors comprising a so-called Bianchi 
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field are given in tensorial notation by 

&ij = ^ioj) Aj = -^VihkVjimR hklm , H,ij = -N 1 r) ihk RQj, Bji = -N 1 r] ihk R 1 ^. (4.7) 

For any Robertson- Walker spacetime, ds 2 = —dt 2 + a 2 (t)da 2 , the norms of the magnetic 
parts, \H\, \B\, are identically zero while \E\ and \D\, the norms of the electric parts, 
reduce to the forms 

\E\ 2 = 3(a/a) 2 and \D\ 2 = 3 ((a/a) 2 + k/a 2 ) 2 . (4.8) 

Therefore the Bel-Robinson energy becomes 

B(t) = ^ (\E\ 2 + \D\ 2 ) , (4.9) 

where C is the constant volume of (or in in the case of a non-compact space) the 3- 
dimensional slice at time t. Thus we find that 

B(t) ~ k 2 u (t) + kl(t), (4.10) 

where k u , k a are the principal sectional curvatures of this space. Using this result, it is 
not difficult to show the following theorem concerning the completeness of Friedmann 
cosmologies following directly from bounds on their Bel-Robinson energy, and an asso- 
ciated minimum radius. 

Theorem 4.4 A spatially closed, expanding at time t*, Friedmann universe that satisfies 
7 < Bit) < r , where 7, Y are constants, is causally geodesically complete. Further, there 
is a minimum radius, a min > A -1 / 2 , A depending on T, and these universes are eternally 
accelerating (a > 0). 

The ensuing classification meeting all the different criteria discussed previously results 
from the possible combinations of the three main functions in the problem, namely, 
the scale factor a, the Hubble expansion rate H and the Bel Robinson energy B. The 
singularity types will by necessity entail a possible blow up in the functions \E\, \D\. If 
we suppose that the model has a finite time singularity at time t — t s , then the possible 
behaviours of the functions in the triplet (H, a, (\E\,\D\)) are as follows: 
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51 H non-integrable on [ti,t] for every t > ti 

52 H — > 00 at t s > ti 

5 3 H otherwise pathological 

N 2 a — > a s 7^ 

£>i l^l — > 00, |D| — > 00 
Z? 2 l-^l < 00, |.D| — > 00 
B 3 \E\ — >■ 00, < 00 
£4 |£| < 00, |.D| < 00. 

The nature of a prescribed singularity is thus described completely by specifying the 
components in a triplet of the form 

{S i ,N j ,B l )i 

with the indices I taking their respective values as above. Here category S monitors 
the asymptotic behaviour of the expansion rate, closely related to the extrinsic curvature 
of the spatial slices, iV that of the scale factor, describing in a sense what the whole of 
space eventually does, while B describes how the matter fields contribute to the evolution 
of the geometry on approach to the singularity. We know from the completeness theorems 
that all these quantities need to be uniformly bounded to produce geodesically complete 
universes. Outside complete universes, the whole situation can be very complicated and 
the classification above exploits what can happen in such a case when we consider the 
relatively simple geometry of isotropic cosmologies. 
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